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STATUS REPORT 


Task $41 is composed of two parts: (l) analysis and design studies related to 
the NAS Extended Operating Configuration (EOC), and (2) Computational 
Chemistry. 


During the first half of 1987, Dr. Levin served as a member of an advanced 
system planning team to establish the requirements, goals and principal 
technical characteristics of the NAS EOC. As a result of earlier studies in this 
area, he presented a paper entitled, “Scaling of Data Communications for an 
Advanced Supercomputer Network,” at the Third International Conference on 
Data Communications Systems and their Performance, in Rio de Janeiro (June 
22-25). The paper was published in the Proceedings of that conference. 


A major milestone was reached in a study being conducted by Dr. Levin of 
RIACS, together with Dr. Partridge and Dr. Stallcop of NASA/Ames. The high 
temperature transport properties (such as, viscosity, thermal conductivity, etc.) 
of the major constituents of air (Oxygen, Nitrogen) were correctly determined 
for the first time. The results of prior ab initio computer solutions of the 
Schrdoedinger equation were combined with the best available experimental 
data, to obtain complete interaction potentials for both neutral and ion-atom 
collision partners. These potentials were then used in a computer program 
developed by Dr. Levin to evaluate the collision cross-sections from which the 
transport properties could be determined. The results are contained in a paper 
entitled, “High Temperature Transport Properties of Air,” by Levin, Partridge 
and Stallcop, which was presented by Levin at the Twenty-Second AIAA 
Thermophysics Conference in Honolulu (June 8-10). These results are of 
importance for the Advanced Orbital Test Vehicle (AOTV) and the Bational 
AeroSpace Plane (NASP). 



PAPERS PUBLISHED 


Eugene Levin, “Scaling of Data Communications for an Advanced Supercom- 
puter Network,” Proceedings of the Third International Conference on Data 
Communications Systems and Their Performance, June 22-25, 1987. 

Levin, Partridge and Stallcop, “High Temperature Transport Properties of Air,” 
presented at the Twenty-Second AIAA Thermophysics Conference, June 8-10, 
1987. 
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The goal of NASA’s Numerical Aerodynamic Simulation (NAS) Program is to 
provide a powerful computational environment for advanced research and 
development in aeronautics and related disciplines . The present NAS system 
consists of a Cray 2 supercomputer connected by a data network to a large 
mass storage system , to sophisticated, local graphics workstations and by remote 
communications to researchers throughout the United States . The program 
plan is to continue acquiring the most powerful supercomputers as they become 
available . This paper describes the implications of a projected 20-fold increase 
in processing power on the data communications requirements. 

1. INTRODUCTION 

The Numerical Aerodynamic Simulation (NAS) Program was initiated by NASA to 
establish a national resource for advanced research and development in aeronautics and 
related disciplines. To achieve this goal the NAS Program is to, "act as the pathfinder 
in advanced, large-scale computer system capability...” [1,2.3, 4]. The first major mile- 
stone has been achieved and the initial configuration (Figure 1) is now operational at 
the NASA Ames Research Center. The centerpiece is a Cray 2 supercomputer with 256 
million (64-bit) words of memory and a sustained performance rated at 250 MFLOPS 
(250 Million FLoating-point OPerations per Second) as measured on benchmark tests 
for optimized large-scale computational aerodynamic codes. 

The need for much more powerful processors can be seen from Figure 2 [l, 2, 3] which 
depicts the estimated speed and memory requirements for various approximations to the 
governing equations and three levels of geometric complexity; an airfoil, a wing and a 
complete aircraft. Note, for example, that if viscosity effects are included by using the 
Reynolds-averaged Navier-Stokes equations, a three dimensional solution for a wing 
requires about 100 times the computing speed of a comparable inviscid solution and 
only now has it become feasible to perform highly repetitive design optimization studies 
for such cases. Furthermore, if still more realistic large eddy effects are considered, a 
further factor of about 1000 is required for runs of the same duration. Finally, Figure 2 
indicates that (with 1985 algorithms), a single 15 minute run including large eddy 
effects for a complete aircraft would require computer speed in excess of 10 12 floating 
point operations per second and a random access memory of about 10 11 bytes! 

* Work supported in part by Cooperative Agreement NCC 2-387 from National 
Aeronautics and Space Administration to Universities Space Research Association. 
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Figure 1. Initial Operating Configuration 
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Figure 2. CPU Speed and Memory Requirements 



In consideration of these computational needs, the NAS Program plan is to continue to 
acquire the most powerful supercomputers as they become available. By 1988 it is anti- 
cipated that a "one GigaFLOPS" computer with 4 times the .si/*/ n: uv.d speed of;. Gray 2 
will be obtained and by 1990 an additional supercomputer with 1G times the sustained 
Cray 2 speed (4 GigaFLOPS, i.o., 4000 MFLOPS) will be added. To assure that this 
increased power can be fully utilized, it is essential to examine the total supporting 
system-level infrastructure. In particular, it is critical to provide sufficient capacity for 
the very large data files characteristic of fluid dynamics computations and to scale the 
bandwidth of the communication systems to handle the increased traffic. 


2. RELATED WORK 

Several authors have recently considered the "balances" needed between computer 
speed, storage requirements and data communications. Kung [5] presents a model of 
balanced computer architectures for particular classes of computation; however his work 
is primarily concerned with the characteristics of t lie processor itself and not the total 
environment. The same limitation is true of "Amdahl’s rules of thumb" (see [ 1 3 j) which 
Worlton [G] states in the form, "One byte of main memory is required to support eacli 
instruction per second," and, "One bit of I/O is required to support each instruction per 
second." 

Thorndyke [7] docs consider the support environment for the ETA-10 supercomputer 
and views the mass storage subsystem as part of a memory hierarchy consisting of cen- 
tral processor memory, shared memory, local disks, and mass storage. He concludes 
that the capacity ratios between each level of this hierarchy should increase by a factor 
of 1 G: 1 . He also proposes that data communication rates be matched to disk transfer 
rates of 10 Megabytes/sec. Ewald and Worlton [8] note that the Cray XMP/48 also 
exhibits a 1G:1 ratio between local "disk" storage (Solid State Disk, in this case) and 
main memory. Furthermore, they note that historically the requirements at all levels of 
the storage hierarchy have been roughly proportional to the speed of the computer. For 
future scaling, they propose that on-line disk capacity requirements should grow at 
about 2/3 of the performance growth of the supercomputer and that transfer rates be 
increased to achieve a balanced system. 

The 198G work of Wallgren [9] uses similar theoretical scaling laws based on past and 
current experience to project the supporting environment needed for future supercom- 
puters. His primary focus is on storage and not on data communications requirements. 
He correctly notes that the results are dependent on the assumed system architecture 
and on the usage profile,. Wallgren's extrapolations extend over more than a decade in 
time and a factor of as much as 10* increase in supercomputer speed. 


3. SCOPE 

The present, study examines the impact on the data communications of a 20-fold 
increase in computer speed over a 3-4 year time period and specifically assumes an archi- 
tecture (Figure 3) which is structurally the same as the existing initial operating confi- 
guration of Figure 1. Furthermore, the user population is well defined since the NAS 
Program is not a general purpose scientific computing center but is devoted specifically 
to aeronautics studies and applications; the program plan calls for some 90% of the 
total time to be used for computational fluid dynamics and aerodynamics. This permit- 
ted a detailed workload model to be developed. 
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Figure 3, 1990 Model Configuration 


4. ANALYSIS OF DATA COMMUNICATIONS 

The 1990 system architecture shown in Figure 3 consists of a primary high speed data 
network between two supercomputers and a central mass storage system; a secondary 
network for communication between local workstations and the remaining subsystems; 
and a remote communication subsystem. The primary focus of the study was to deter- 
mine the requirements for the "backbone” high speed data network. for the movement of 
large files (typically from 5-80 million words) to and from the high speed processors. 
The remote communication subsystem and the local area network serving the worksta- 
tions are also scheduled for upgrade, however their sizing was determined by factors 
other than the speed of the supercomputer(s) . 

4.1. REMOTE COMMUNICATIONS 

The NAS remote communications system currently supports Arpanet/Milnct, NFSnet, 
NASnet, and the NASA-wide Program Support Communication Network (PSCN). 
Access to the NAS network from remote sites is shown in Figure 4. Vitalink bridges 
manage the inter-network connection and monitor all traffic on the ethernct local area 
network passing only those messages that have remote destination addresses to the 
appropriate Vitalink unit. 

Doth terrestrial links (at 5G and 244 Kilobits/sec) as well as satellite Ti links (at 1.544 
megabits/ sec to the NASA centers) arc provided. At present, over 20 remote sites (with 
approximately 100 users) have been activated and additional sites will be added to 
expand the user community. 
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Figure 4. NAS Remote Communications 







The long term goal is to provide equivalent services to both local and remote NAS users. 
This objective is constrained by technology and funding limitations. T^)grade to T2 
rates (G.2 megabits/ sec) for selected NASA centers is planned during the 1938/1989 
time period. Additional improvements under investigation include implementing class 
of service protocols (such as distinguishing bulk file transfers from interactive traffic) 
and techniques for reducing the volume of remote data communications (such as remote 
block editors and distributed graphics services). 

4.2 LOCAL WORKSTATION COMMUNICATIONS 

The principal local users access the NAS system through powerful graphics workstations 
(Silicon Graphics 2500 Turbos). The workstat ion is an essential element of the system 
to permit graphical analysts and interpretation of the extremely large output files gen- 
erated by computational fluid dynamics batch programs. It also enables the user to 
work interactively with the supercomputer. For example, a user may designate loca- 
tions in the vicinity of an aerodynamic surface on his display, the flow patterns 
representing particle trace's from those locations are computed on the supercomputer 
and the results returned to the user's display. The data communications services to the 
workstation must be sufficient to enable the user to make maximum use of both the 
batch and interactive workstat ion/suporcoinputer capabilities. 

There are three modes in which graphics can be displayed at the workstation. It is pos- 
sible to send down the solution files (or large subsets of these files) for computation and 
display at the local workstation. Although the workstation lias the processing power 
comparable to a \ ax 11/780 and the special purpose graphics hardware cxj>cditcs the 
generation and transformation of vectors and polygons, there arc many limitations to 
this mode of operation including the I/O bandwidth of the workstation and limitations 
of main memory and disk space. The preferred modes use the supercomputer for the 
computationally intensive tasks and send down either display lists ox pixel data to the 
workstation. Both modes may be used, however with the current workstations it is diffi- 
cult to use the pixel mode of interaction effectively. Future workstation upgrades will 
alleviate this limitation and the 1990 model assumes that both the display list mode and 
the pixel mode of interaction will be used extensively. 

Based on the above description of operational use, the data communications require- 
ments to the workstation are primarily sized according to the graphics needs of the user 
and the projected capabilities of future workstations. The present effective rate of 
approximately 2 Megabits per second is appropriate for the current workstation receiv- 
ing display lists (typically 2.G Megabits in size) from the supercomputer since the 
transmission time (1.3 seconds at 2 Megabits per second) is only a fraction of the time 
needed for the workstation to generate the display (approximately 5 seconds). Some 
representative display capabilities at. 8 Megabits per second (the minimum planned for 
1990) arc shown in Table 1. For example, a minimal pixel file of 8 Megabits (display 
resolution of 1,000 x 1,000 with 8 bits/pixcl for simple color and shading) ran be 
transmitted in 1 second and the time for the workstation to generate a display from 
such a file is negligible. A very high quality pixel file of 9G Megabits would require 12 
seconds of transmission time but this would typically be used for presentation quality 
graphics only. (For comparison, the data rate needed to support real time, interactive, 
color graphics with animation at 30 frames/sec. would he about 800 megabits/sec. The 
NAS Program is funding a prototype of such a bus and considering prototyping an 
advanced graphics workstation with a massive frame buffer and very high resolution, 
however this is not included in the 1990 system model.) 



Table 1. Display of Color Images 


Type of Output 
from Supercomputer 
to Workstation 

Typical 
File Size 

Channel Time 
@ 8 Mbps 

Time for 
Workstation to 
Generate Display 

(8000 Polygons, 

Current Wo 
includes Hidden 5 

rkstation 
Jurface Remova 

1, Flat Shading) 

Solution File 

400 Megabits 
(near storage 
limit) 

30 seconds 

= 25 seconds 

Display List 

2.6 Megabits 

0.3 seconds 

~ 5 seconds 
( ~ 1 second if no 
hidden surface 
removal) 

Future Workstation 

Minimal Pixel File 

8 Megabits 
(lk x lk x 8 bits) 

1 second 

« 1 second 

Fully Rendered, 
Very High Quality, 
Pixel File 

96 Megabits 
(2k x 2k pixels 
24 bits/pixel) 

12 seconds 

« 1 second 


4.3. HIGH SrEED "BACKBONE" COMMUNICATIONS 

The high speed data communicat ion requirements were derived by analysis of a model of 
t lie expected workload and checked using a discrete simulation program. The workload 
model assumed that the processing power of 5 Gigaflops available from the two 1990 
high speed supercomputers (Figure 3) was fully utilized. A detailed profile (for compu- 
tational fluid dynamics) of types of user tasks and expected frequencies had been 
developed over a six year time period and was updated in 198G to reflect current algo- 
rithms, projected increased interactive usage and substantially increased use of graphics 
[ 10 , 11 , 12 ]. 

The model contained over two hundred parameters including: numbers and types of 
local and remote users, number of host processors, protocol delays, amount of disk 
storage attached to the supercomputers, distribution of batch and interactive work, fre- 
quency of task execution, probability of abortive runs, etc. Scripts were developed to 
represent characteristic delays associated with "think time" to separate user-initiated 
sequential processes. These asynchronous processes compete for system resources. The 
high speed data communications capability was initially taken as unbounded and the 
workload was progressively increased until the full capability of the supercomputers was 
saturated. 

A simplified listing of the classes of work initiated by the users is summarized in Table 
2. The principal execution runs represent various types of numerical aerodynamic simu- 
lations used to solve fluid flow problems. These generally involved repeated iterations of 



















difference equations over a three-dimensional grid which was assumed to consist of one 
million grid points. The simple steady state design simulations used an in vise id poten- 
tial and required approximately 20,000 calculations per grid point of result File. The 
more complex steady state design simulations used the Reynolds- Averaged form of the 
Navier-Stokcs equations and required approximately 000,000 calculations per grid point. 
The comparable unsteady solutions typically required over 4 million computations per 
grid point. The "model day" included 73 hour-long runs of this computation, each 
requiring 1.0 (JigaFLOPS of sustained processing power. From Table 2 it is evident 
that this represented the dominant load on the supercomputers. 


Table 2. Estimated Processor Load of User Initiated Runs 
TASK Q flop/pay 


CODE & PARAMETER PREPARATION 6 6 

PATCH GENERATION ^ 43 

GRID GENERATION 96 

METHOD AND CODE DEVELOPMENT 13,003 

SIMPLE STEADY STATE DESIGN SIMULATIONS 5*732 

COMPLEX STEADY STATE DESIGN SIMULATIONS 31,834 

COMPLEX UNSTEADY STATE DESIGN SIMULATIONS 25 0,978 

FLUID RESEARCH LARGE EDDY SIMULATIONS 53,472 

RESULT EDITING & VIEWING 17,465 

DOCUMENT PREPARATION & USER COMMUNICATIONS * 31 
X 0 T A L 372,719 


In order to relate the computational loads to network traffic, a Workload Data Base 
was constructed which contained detailed estimates of the sizes of the various 
in put /output, files for each fash and a model of how they would be utilized. This pro- 
vided a means for determining the specific files generated and their movement across the 
network. For example. Table 3 shows the network traffic load for the dominant 
Unsteady Design Simulation task. 


Table 3. Network Load for Unsteady Design Simulations 

MB/day 


Sources 
D e s t i n a t i ons^s^ 

MASS 

STORAGE 

HIGH SPEED 
PROCESSORS 

TOTAL 

MASS 

STORAGE 

- 

4,428 

4,428 

LONG-HAUL 

COMMUNICATIONS 

0.2 

2,048 

2,048 

WORK- 

STATIONS 

0.1 

1,810 

1,810 

HIGH SPEED 
PROCESSORS 

1,555 

- 

1,555 

TOTAL 

1,555 

8,286 

9,841 
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The source of most of the data communications load comes from the output of large 
Tiles from the the high speed processors to their local disks and the subsequent move- 
ment of these files between the supercomputers and the central mass storage system. 
Figure 5 shows the results of the model study for this traffic. There is a predicted net 
accumulation on the local disks of 42 Gigabytes per day which must be moved (by 
system-managed file migration) to the central mass storage system so that the super- 
computers 1 disk space will be available for current work. The two supercomputers of the 
1990 system were modelled with a total of 400 Gigabytes of local storage. Since more 
than half the available 400 Gigabytes is needed for current working space each day, the 
average age of files before automatic migration to central storage was found to be about 
3.5 days. This additional load (of automatic migration) on the data communications 
can, however, be distributed over the periods of lightest user activity. Typically the 
nighttime workload consists of large batch jobs with little interactive traffic. Under 
these conditions the supercomputers are fully saturated but the network is only lightly 
loaded. -v 

Note: Supercomputer Disks 
Daily Net Accumulation = 42 GB 


Data Created on HSPs 
165 GB/Day 



Figure 5. Daily File Activity 


Figure 6 summarizes the major result of the study and illustrates both the hourly distri- 
bution of traffic as well as the average daily traffic including migration. To avoid signi- 
ficant queueing delays during times of peak activity, and to handle periodic bursts 
within an hour efficiently, the bandwidth for the 1990 system should be sized well in 
excess of the peak rates shown on Figure 6. The design value recommended was 100 
megabits per second [14]. This was regarded as a reasonable requirement in view of pro- 
jected technological advances. 
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5. CRITIQUE OF RESULTS 

An essential (but often neglected) step in studies of this type is to estimate the sensi- 
ti\ ity of t lie results to t lie assumptions inherent in the model. The succeeding sections 
provide estimates of the effects of some of the major assumptions on the projected data 
communications requirement. 

5.1. USE OF GRAPHICS 

The need for sophisticated graphics to analyze the very large data files typical of com- 
putational fluid dynamics is unquestioned and a significant allowance for this type of 
processing was included in the model. However, this was a projection and might not 
conespond to the actual future usage. Among other factors it is dependent on the capa- 
bilities of future workstations. 

The model provided a breakdown into various types of output. The three principal 
classes created were raw solution results, graphics display lists and pixel files. It was 
found that a very useful measure for scaling data communications with processing 
power was the ratio of bytes of output to MFLOPs (millions of floating point opera- 
tions) required to generate that output. This characteristic output ratio will be desig- 
nated /? M with units of Bytes/MFLOP. The importance of this parameter is that if the 
only deviation from the assumed model is processing power, the results can (within rea- 
sonable limits) be scaled directly. Furthermore, if the mix of classes of output varies, it 
is possible to estimate the effect on this critical parameter and hence (within modest 
bounds) on the final results. 

Over the distribution and frequency of tasks initiated by users, the model showed that 
raw result files (including debug runs) were found to generate /? M = 130 B/M (Bytes per 
MFLOP). In contrast, the generation of graphics display lists generate /? M ^8,000 B/M 
and pixel files generate /? M «5000 B/M. In the model, 5% of the total available process- 
ing power (FLOPs) was devoted to graphics. Approximately 75% of the graphics 




FLOPs were utilized for pixel files, however much more processing is required for pixel 
files. Hence, this distribution corresponded to approximately three graphics display lists 
for each pixel file produced. Since the supercomputer power of the 1990 model system 
was sized at 5 G FLOPS, the 5% utilized for graphics corresponds to 250 MFLOPS, 
equivalent to the total sustained processing capability of the 4 CPU’s of a CRAY 2! 

For the assumed usage distribution, the overall value of the characteristic output ratio 
for the 1990 model was found to be /?^~40() Dyt.es/Mfiop. Figure 7 shows the variation 
in this parameter with changes in the relative amount and type of graphics usage. The 
abscissa represents the proportion of processing power used for graphics processing 
(both display list and pixel data) and the one-parameter family of lines represents vari- 
ous mixes of pixel output- and display list output. The reference point corresponding to 
the model is marked. As an example, the calculation of the design 'point (using the data 
values of the paragraph above) is shown below: 

0.95(130) + 0.05[0. 75(5000) + 0.25(8000)]=411 Bytcs/Mflop. 

It may be noted that increased use of pixel data relative to display list output has a sig- 
nificant effect in reducing fi^. Display list output with a /? M of 8,000 B/M (compared 
with pixel file output with a of 5000) utilizes less supercomputer processing power 
but places a greater load on the data communications (and oil the workstations). 



relative amount of HSP computing power used for graphic processing 

Figure 7. Influence of graphics processing 



5.2. THE USER PROFILE 


The wondoad model was specifically designed for solving problem- in computational 

fluid dynamics. Table 4 summarizes a 198G study of major users of supercomputers at 
the NASA Ames Research Center. The first two columns represent supercomputers of 
the Ames Central Computing Facility and only the Cray 2 is a NAS supercomputer. 
Furthermore, the sampling period shown in Table 4 occurred shortly after the Cray 2 
was initially installed, in a pro-operational time period, and does not represent the 
inf ended operational distribution of users. Nevertheless, it appeared that computational 
chemistry utilized large amounts of supercomputer time and might represent a future 
NAS system user that was not reflected in the workload model. 


Table 4. Ames Vector Processor Utilization 
Percentage of User Time by Organization 


USERS 

CYBER 
205 
FY 85 

CRAY 
X-MP 48 
3/86-10/86 

CRAY 

2 

3/86-10/86 

Remote 


1 

i 

i 

i 

i 

a 

a 

a 

36.8 % 

Advanced 

Aerodynamics Concepts 


7.1 % 


Applied 

Computational Fluids 

0.1 % 

24.1 % 

7.0 % 

Experimental 

Fluid Dynamics 

0.7 % 

6.1 % 


Computational 
Fluid Dynamics 

25.0 % 

24.1 % 

13.8 % 

Computational Chemistry 

67.9 % 

16.0 % 

19.2 % 

1 

TOTAL 

93.7 % 

77.4 % 

76.8 % 


Much of the current work in computational chemistry at NASA Ames Research Center 
is devoted to determining the electron density distribution of molecular structures by 
solving the Schroedinger equation. The methods used involve evaluation of many 
multi-dimensional integrals and manipulation of very large matrices to find eigenvectors 
and eigenvalues. These methods are quite different than those used in computational 
aerodynamics studies. Hence the output of such runs might produce larger or smaller 
files than those used in the model. 

Consequently, a study was made of the character of the output files generated by these 
supercomputer users. The results are shown in Figure 8 in terms of the parameter 
The prediction of the 1990 model study is shown for comparison. The characteristic 
output of computational chemistry was found to be comparable to the major computa- 
tional aerodynamics users. Hence, even if computational chemistry were in the future to 
become a significant component in the mix of NAS users, it was concluded that this 
would not invalidate the results. 
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Figure 8. Comparison of Current Usage 
and 1990 Model 


It is also of considerable interest, to note the relationship of the predicted 1990 usage 
compared to actual 198G usage. It is not surprising that the predicted value of /? M 
(~400 Bytcs/MFLOP) is almost twice the 1986 measured values in Figure 8 since exten- 
sive use of graphics was just beginning. Figure 7 confirms that increased graphics use 
results in an increase in /? M . 

5.3. AMOUNT OF LOCAL DISK STORAGE 

Locality of files increases the hit rate and reduces the volume of data traffic (sec [15]). 
If local storage were unlimited, there would he no need to move files to a central storage 
facility. If files could be maintained local to the supercomputers for about 10 days (two 
working weeks) instead of 3.5 days, it was estimated that the recall load would he 
reduced from 20 Gigabytes per day to less than 10 Gigabytes per day. Aging would 
also increase the number of files deleted and decrease the amount of manual migration 
shown in Figure 5. All of these effects decrease the traffic on the network. It was 
estimated that increasing the local disk storage by an additional 250 Gigabytes could 
reduce the traffic between the supercomputers and central storage by 30-50%. 

5.4. OTHER CONSIDERATIONS 

The model assumed a typical grid size of 10° points. Increasing the geometrical refine- 
ment by using more grid points does not have any affect on /? M provided the computa- 
tions needed per grid point remain constant. However, in many cases not all the final 
data needs to be saved. The finer grid not only yields results at more spatial locations 
but also permits a more accurate approximation to the differential equation. It is often 
only necessary to save the results on a much coarser grid (except for critical regions of 
interest) and recompute the interior points if needed. For example, saving every fourth 
grid point in each dimension would reduce the output by a factor of 64. Similarly, it 
may not always be necessary to save 64-bits of precision in the output data. 



The realism of the physics was discussed briefly in connection with Figure 2 where it 
was noted that solving the more complex equations such as large eddy simulation 
involved an enormous increase in the number of compulations. There would certainly 
not be a proportionate increase in the output, hence the trend towards removing the 
approximations would result in a sharp decrease in the ratio /? M and alleviate the data 
communications load. This was a factor in the 1990 model and will become increasingly 
significant as more powerful supercomputers become available. 

It was assumed that, the algorithms used in the 1990 time frame would be the same as 
those presently known or in use. This is quite unlikely since algorithm development is a 
fertile, on-going activity. Prior trends have been in the direction of progressively 
decreasing /? M , but it is impossible to predict the implications of undiscovered algo- 
rithms. If there is a radical change from deterministic approaches to statistical tech- 
niques (for example, gas lattice automata methods, [ 1 G] ) this would clearly require much 
more computational effort without a corresponding increase in output, resulting in a 
sharp decrease in 0 M . 


G. SUMMARY 

The results of an analysis of the 1990 ’’backbone” data communications requirements 
for the advanced Numerical Aerodynamic Simulation Program showed that the expected 
workload required an effective bandwidth of 100 megabits per second. 

A critical examination of the various assumptions inherent in the model indicated that 
100 Mbps was a safe, conservative result. 

Among the most sensitive assumptions was the projected amount and type of interac- 
tive graphics expected to be used. 

Increasing the amount of disk storage local to the supercomputers would result in a sig- 
nificant decrease in the data communications requirements. 

The characteristic output ratio (representing a measure of bytes of output relative to 
megaFLOPs of processing) made it possible to estimate the effects of (modest) varia- 
tions in the model (assumptions on the data communications requirements. 
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ABSTRACT 

We have developed a general computer code to allow calculation of atom-atom and 
ion-atom transport collision integrals from accurate potential energy curves described 
by a set of discrete data points (e.g., RKR turning points and energies derived from 
molecular structure computations) for a broad range of scattering conditions. This 
code is based upon semiclassical approximations that properly account for quantum 
mechanical behavior such as tunneling effects near a barrier maximum, resonance 
charge exchange, and nuclear symmetry effects. We have determined transport colli- 
sion integrals for N— N, O-O, N + — N, and 0 + — O interactions from complete sets of 
accurate potential functions derived from combined experimental and ab initio struc- 
ture calculations. For the 0 — 0 case, this includes results for excited states; i.e., all 
possible collisions involving 0( 3 P), 0( 1 D), and 0( ] S). Our calculated values of the 
N^— N and 0^ — 0 resonance charge exchange cross section Q ex agree well with meas- 
urements from beam experiment that are available at high energies where the diffu- 
sion cross section Q d satisfies Q d ~ 2Q ex . Tables of collision integrals and plots of 
viscosity, thermal conductivity, and diffusion coefficients that are needed to determine 
the transport properties of high temperature air are presented. 
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I. INTRODUCTION 


Transport properties involving high temperature air are of increasing importance. In 
particular, they are needed for a real gas analysis of the nonequilibrium shock layer 
flows about Aeroassisted Orbital Transfer Vehicles (AOTV’s) during their aerobrak- 
ing maneuvers at high altitudes in the earth’s atmosphere and for design studies of 
the National Aerospace Plane (NASP). 

The determination of transport properties requires knowledge of the potential energy 
curves for the various interacting species. Existing calculations have, by and large, 
been based on incomplete potential information; i.e., generally only the RKR 
(Rydberg-Klein-Rees) curves near the well minimum for a few bound states and often 
meager potential data on the higher lying molecular states. Furthermore, the method 
of calculation of the transport properties was based upon fitting the potential data by 
adjusting the parameters of an assumed analytical form (e.g., Morse potential). Such 
analytical forms fail to exhibit the correct long-range interaction potential. These 
various omissions and approximations limit the reliability of previously computed 
transport properties. 

The present work incorporates three major improvements: accurate representation of 
the interaction curves for a very large set of interacting atomic and ionic species, use 
of spline curves to fit the available potential data with no a prion constraints as to 
the functional form, and incorporation of certain quantum effects. We have carefully 
combined the available experimental data with interaction energies obtained from ab 
initio structure calculations to construct the best representation of all the potential 
curves for O — O, N-N, N + -~N, and 0 + -0 when the collision partners are in their 
ground states. In addition, we have included the excited states of 0-0 to allow the 
assessment of their contribution to transport properties under nonequilibrium condi- 
tions. In all cases the long-range interactions were also correctly represented to cover 
the entire range of separations necessary for the accurate calculation of transport pro- 
perties. 

The potential data was constructed from complete sets of accurate ab initio calcula- 
tions. For example, our computed potential curves for N + — N, and 0 + — O interactions 
[l] agree to within 0.2 eV with the RKR curves in the region of the potential 
minimum. Our calculated long range potential energies have been applied to deter- 
mine resonance charge exchange cross sections [2] that agree w r ell with the 
N + — N, and 0 + — O cross sections measured by beam experiments. In the case of 
N + — N interactions we have found that potentials which were previously thought to be 
entirely repulsive, were instead, in fact, bound [3], and hence could affect the predicted 
transport properties. Our calculation of complete sets of accurate ion-atom potential 
energy curves for the N + — N, and 0 + — O interactions permits the determination of 
reliable transport properties for these species for the first time. 

A second feature of the present work is the use of our new computer code, "SKATR" 
that enables the calculation of collision integrals for a broad range of scattering condi- 
tions from the actual potential energy as specified by a discrete set of data points, 
rather than approximate analytical fits. The program is structured to fit a spline 
curve through the discrete data and subsequently to join smoothly with a correct 
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specification of the long-range interaction. The program also incorporates the quan- 
tum mechanical effects of the semiclassical formulation described below. 

We have replaced the non-physical singularity of the classical deflection function at 
the potential energy barrier maximum by a uniform phase shift approximation due to 
Stallcop [4], which includes the quantum effects of barrier tunneling and resonance 
scattering arising from the metastable energy levels of the inner well. Other quantum 
mechanical behaviour taken into account by our approach includes resonance charge 
exchange, which dominates the like ion-atom diffusion process at high energies, and 
nuclear symmetry effects which are significant at low energies. Application of the 
semiclassical approximation should yield accurate cross sections even at low collision 
energies. This prepares a foundation for future studies involving very light collision 
partners such as hydrogen (important for the propulsion system of the NASP). 
Furthermore, the semiclassical approach allows us to calibrate the accuracy of the 
classical method (which requires much less computation) and to ascertain its range of 
applicability. 

The mathematical formulation including the semiclassical approximation is described 
in Section II. In Section III, we present the details of the construction of the accurate 
interaction curves used as input. A brief description of the SKATR program is 
presented in Section IV. The calculated results are shown and discussed in Section V, 
followed by concluding remarks and observations in Section VI. 


II. METHOD 

According to Chapman-Enskog theory [5], the transport properties of a dilute gas can 
be expressed in terms of the transport collision integrals, 


^n.s ( T ) = F ( n > s ) f e *V S+3 Qn (kTx 2) dx, 


( 1 ) 


where F(n,s) = 


4(n+l) 


Jr(s+1)! [2n+l — (— l) n ] ’ 


k is the Boltzmann constant, T is the temperature, o is the hard sphere collision diam- 
eter, and <r 2 n is the mean collision integral. 

For a gas mixture to first_order,_or a pure gas to second order, the transport collision 
integrals required are: a 2 Cl 1 l5 a 2 Oj 2 , a 2 Q 1 3 , <T 2 n x 4 , o 2 Vl l 5 , «x 2 0 2 2> 2 3> 4> an< i 

<r 2 fi 3j3 [6]. 

For example, for a pure gas to first order, 

Viscosity: i] — 2.6693xl0 -5 — , g cm -1 sec 

a 1^2, 2(T) 
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Thermal Conductivity: 


A = 0.19891 


(T/M)* 
<^2,2(T) ’ 


meal cm 1 sec 1 K 1 


Self— diffusion Coefficient: 


D = 2.628x10 


— 3 (tVm)* ^ 

p^rii ^T) ’ 


cm 2 sec 1 


where M is the molecular weight in g/mole, a 2 Q is in square Angstroms, T is the Kel- 
vin temperature, and p is the pressure in atmospheres. 


For a gas mixture, the corresponding diffusion coefficient becomes 

Binary Diffusion Coefficient: D:: = 2.628xl0~ 3 — , cm 2 sec -1 

P a 2 n M (T) 


where fi = 


MjMj 

Mi + Mj * 


The pure gas viscosities, thermal conductivities, and diffusion coefficients can be used 
together with the binary diffusion coefficient and species concentrations to determine 
the viscosity and thermal conductivity of a gas mixture, such as air. These resultant 
transport properties may then be used as inputs to vehicle flow field calculations. 


The evaluation of the collision integrals of Equation (l) requires knowledge of the col- 
lision cross sections Q n (E) for a broad range of values of the parameter, E. These 
cross sections can be calculated from classical scattering theory or determined from 
quantum mechanical scattering phase shifts. At high energies, the classical scattering 
approximation is valid and the cross sections may be evaluated from 


00 

Q n (E) =2 7T f 6(1 — cos"x) d6, where 

J o 


( 2 ) 



dr 

7¥i7) 


In these equations, b is the impact parameter and x * s the angle of deflection (See Fig- 
ure l); r t is the classical M turning point”, the largest root of 


F( r ) = 


E r 2 


% 


= 0, 


and V(r) is the interaction energy of the collision partners. 
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t : 




scattering center 


Figure 1. Scattering Geometry 


The corresponding quantum theory expressions for the transport cross sections, Q r 
for ion-atom collisions are shown below [7]: 


1 + , r\ i ( 


2s -r A 


Qn 2 [Qn + Qn 1 + 2(2s + jj 


[Qn + - Qn']. 


(0 


where A = 0, 1, 2, ....for the £, II, A, ...etc., states, 

{ 1 for nitrogen, 

n C 

0 for oxygen. 

For odd values of n, 


Qn + = 7T S f Qodd(g. u ^) + Qeven( u .g^)l 

* Vm>0 L J 


(5a) 


Qn = E f Qevenls."^) + QoddKs; 1 ') 1 

k •* L J 


(5b) 


*^>0 



- 5 - 


and, for even values of n, 


Q„ + = 


4tv 


j QevenfS^Si^ 7 ) 4r Qoddfaj^i^)! 

^v«>0 L J 


(6a) 


Qn _ = 


4 7T 


S f Qodd(S’Sj^) *b Qeven(^’^’^) 1 

L J 


^«v eO >0 


(6b) 


where, 


Qeven(f»’ u > 1/ ) E a n is s * n (*?/+ u Wl ) 


QoddCe. 11 ;^) = E a n«/ - *?“) 

/ 0 dd 


In these expressions, g and u refer to the gerade and ungerade states respectively, and 
/ is the angular momentum quantum number. 

The coefficients a n ^ are determined by recursion from 


(21 + 1) x" P,(x) = £ ai P, + „(x) 


where P;(x) are the Legendre Polynomials [8]. The charge exchange cross section is 
given by 


r\ _ ^2, (21 + 1) . 2/ e u\ 

Qex = TT E Lsin (W ~ Vl ) 

k 1=Q 4 


(7) 


The transport cross sections for neutral atom-atom collisions were calculated from 


" TT S E a^sin 2 (»?,+*,— »?/), 

*■ »/>0 1=0 


( 8 ) 


even 


where the allowed values of v are < > according to the even or odd parity of n. 
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At high energies, the scattering phase shifts rji are approximated by the JWKB for- 
mula. 






dr 


The lower limit r 3 of the first integral is the largest root of 


G(R) = 1 


v(R) 

E 



= 0 . 


( 9 ) 


In Equation 9, /?=(/ + 2 / 2 ) /k can be identified with the classical impact parameter b. 

At low energies, we apply a semiclassical method that accounts for tunnelling through 
the barrier of the effective potential energy (see Fig. 2). The uniform phase shift 
approximation developed by Stallcop [4] also accounts for resonance scattering associ- 
ated with metastable energy levels of the inner potential well and virtual energy levels 
above the barrier minimum. 



Figure 2. Effective Potential Energy, V e 
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For energy levels below the barrier, 
f] = *7j w ](b + P + arctan 


{ tan( * + ■ where 


0.5 arg rl ■/, + —] 

[-fin 

I 

1 

| 2ir 

G7T 


r=0.5 1n{ r (‘l, + i-)lV2i\ - In | | , 

7 r I 27r en 


$ = k j [G(r)]^ dr, 


Tl 


I = k f | G(r) | % dr , 

rj 

and r 1 r 2) r 3 are the ordered roots of G(r)=0. 

For energy levels above the barrier, there is only one root, r l5 and r) takes the form 

sin 2(<f> + p) 


V = ^jwkb + 2 p -arctan 

I'm 

^ = k / [G(r)] % dr, 


vTT e 21 + cos 2{<f> + p) 


}■ 


where now, 


nr J 

/ il 1 

I , 

i 

p — 0.5 arg T< 


In 


I 

i M 

27T 

e7r 


, and 


I = 


— 7rk 


I 


G 


v/2 \ \/| G"| I r = r „ 


The quantity r=r m specifies the position of the barrier maximum. The expression for I 
is based on a local parabolic expansion near the barrier maximum. 
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III. SOURCES OF INTERACTION ENERGY DATA 

The interaction energies (which correspond to the potential energy curves of the vari- 
ous states of the molecules formed from the colliding atom-atom or ion-atom pairs) 
are needed as input to determine the transport cross sections. We have combined 
experimental data with the results of ab initio molecular structure calculations to 
obtain the best representation of the true interaction energies. 

For example, where available we have used experimental data such as the measured 
dissociation energy [12] and RKR turning points to represent the potential in the 
region of the well minimum. Ab initio molecular structure potential calculations were 
used to extend the experimental data to the adjoining regions by adjusting the ab ini- 
tio results to agree with the experimental data. Since electron exchange interactions 
dominate in this region (i.e., the potential energy behaves roughly like an exponential 
function V 0 e~ aR ), the adjustments were made in such a way as to preserve the value 
of the shielding parameter, a. For small R, where the potential curves are steep, we 
have shifted the ab initio potentials by AR, the difference between the calculated and 
measured equilibrium separation distance R e . Since the measured and calculated 
values of the potential curves are in good agreement (e.g., Ref. [l]), the transport 
cross sections are not very sensitive to the manner of adjusting the ab initio results to 
be compatible with the measured data. 

The ab initio interaction curves have been extended by a proper form of the long 
range potential interaction. For the neutral atom-atom interactions, such as N-N and 
0-0, we use the potential form — C 6 R 6 . The ion-atom N + —N interaction exhibits 
the charge-induced dipole force — l / 2 aTl~* , whereas the 0 + — O interaction behaves like 
a charge-quadrupole potential /^Qm^ 3 • The quadrupole moment is negative for 
£ state interactions, but is positive for II states [13]. The values of the parameters 
C 6 , a and have been determined either from the ab initio results at large R or 
taken directly from other approximations to the Schroedinger equation such as pertur- 
bation expansions [14]. 

We now describe in detail the source of all the potential curves; we begin with N-N. 
The potential wells of the bound states X x £g + and A 3 £ u + of N 2 were obtained from 
RKR data [15,16]. Additional potential information for these states and also for the 
5 £ g + state was derived from the configuration interaction studies of Ermler, McLean 
and Mulliken [17], The repulsive region of the ground state X 1 E g + was extended to 
small R with the help of the generalized valence bond calculations of Dunning, Cart- 
wright, Hunt, Hay and Bobrowicz [18] and Guberman [19]. The 5 £ g + potential was 
extended to large R with the MC-SCF results of Krauss and Neumann [20], The 7 £ u + 
state potential has been computed accurately by Partridge, Langhoff and Bauschlicher 
[21]; these results have been extended to smaller R by adjusting the SCF potential 
determined by Ferrante and Stwalley [22]. The long range calculations of [21] 
strongly support the value of the dispersion coefficient deduced by Zeiss and Meath 
[14]; consequently, we have used their values for all the long range N-N interactions. 
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The sources of the 0-0 interactions will be described next. Saxon and Liu [23] have 
performed configuration interaction calculations for the 62 electronic states of 0 2 that 
dissociate into O atoms in their lowest 3 P, 1 D, and *S states. We find that the values 
of C 6 deduced from their large R results are consistent with those of Neumann and 
Krauss [24]. The Saxon and Liu results have been used directly for collisions involving 
excited state atoms. For ground state 0-0 interactions we have made the following 
modifications to the Saxon and Liu curves. The potential wells of the 
X 3 H g - , a*A g b^g^ - , and the A 3 £ u + states of 0 2 have been determined from RKR 
data [25]. Spectroscopic data for the other two bound states were used to determine 
the parameters of an HH potential function [26,27]. The ab initio results were then 
adjusted to the curve obtained from this potential fit in the region of the well 
minimum in the same manner described above for the RKR curves. 

In addition, we have included corrections for some Rydberg states and shifted energies 
at small R to give better agreement with electron scattering measurements. The Ryd- 
berg character of the 3 II U state at small R has been taken into account using the 
potentials calculated by Buenker and Peyerimhoff [28] and Guberman and Dalgarno 
[29]. Similarly, Rydberg corrections to the 3 IIg and *1^ states have been taken from 
Saxon and Liu [30,31]. Furthermore, the potential curves of the 3 II U , and 3 IT g states 
have been shifted at small R to be consistent with the potentials deduced from elec- 
tron scattering measurements [32,33,34]. 

Our calculations of the N + — N interaction curves using a CASSCF method has been 
described previously [1,2,3]. Primarily the results from the large basis set [2] have 
been used for the transport calculations of this work. We have, however, augmented 
the data from this large basis set with additional data points deduced from the small 
basis set calculations. This w r as accomplished by shifting the energies of the small 
basis set by an amount determined by a parabolic interpolation of the energy differ- 
ences between the two basis sets. These additional data points are especially useful in 
the regions of the potential wells and avoided crossings. 

The wells of the bound states X 2 S g + , A 2 II U and 2 n g were obtained from RKR data 
[25]. The potential curve extracted from spectroscopic data by Hajj [35] was used for 
the B 2 E u + state well region. 

Our calculated potentials for the 0 + — 0 interactions are shown in paper [l] . We have 
modified these curves with RKR data for the bound states X 2 n g A 2 IT Ui and a 4 II u . 
The long range quadrupole parameter Q^j was obtained from the calculated energies 
at large R after first subtracting the single electron exchange contribution [2]. 

A summary listing of all the states for which we calculated transport cross sections is 
shown in Table 1. The weights shown following the state designations represent the 
degree of degeneracy (multiplicity) of that state and must be considered in the deter- 
mination of mean collision cross sections. The potentials employed are contained in a 
report [36] which is available upon request. 
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Table 1. LISTING OF STATES ANALYZED 

0-0 STATES (WEIGHTS) 

GROUND 

EXCITED 

0( 3 P)-0( 3 P) 

0( 3 P)-0( 1 D) 

0( 1 D)-0( 1 D) 

0( 3 P)-0( 1 S) 

0( 1 D)-0( 1 S) 

o( 1 s)-o( 1 s) 

bV (1) 
<X“ (1) 

2 x E g + (1) 
X 3 E g - (3) 
A 3 E u + (3) 
2 3 E u + (3) 

X~ (5) 
l 5 E g + (5) 
2 5 E g + (5) 

IX (2) 

IX (2) 
i 3 n g (6) 
i 3 n u (6) 

X ( 10 ) 

X (io) 

a*A g (2) 
C 3 A u (6) 

X (10) 

B 3 E u - (3) 
3 E g + (3) 
2 3 E u - (3) 
2 3 E g - (3) 
3 3 E u + (3) 
3 3 E g " (3) 

2 3 n g (6) 
2 3 n u (6) 
3 3 n g (6) 
3 3 n u (6) 
4 3 n g (6) 
4 3 n u (6) 
l 3 A g (6) 
2 3 A g (6) 
2 3 A u (6) 
3 3 A g (6) 

X (6) 
X (6) 

*X" (1) 
«X" (1) 

3 J S g + (1) 
4%+ (1) 
(i) 

2 X (2) 
*X (2) 
3 X (2) 

3X (2) 

l 1 A g (2) 
2 1 A g (2) 
3'A U (2) 

X (2) 
X (2) 
X (2) 

3 3 E u - (3) 
4 3 E,- (3) 

s 3 n t (6) 
5 3 n,, (6) 

l^u + (1) 
6 X E g + (1) 
4 x n g (2) 
4 x n u (2) 
2 1 A u (2) 
4 1 A g (2) 

7 x E g - (1) 




N( 4 S°)-N( 4 S°) 

N + ( 3 P)-N( 4 S°) 

0 + ( 4 S°)-0( 3 P) 

*X + (1) 

A 3 E u + (3) 
X + (5) 

X + (7) 

X + , 2 E u + (2) 

4 E g + , X+ (4) 
6 S g + , 6 E u + (6) 
2 n g , 2 n u (4) 
4 n g , 4 n u (s) 
6 n g , 6 n u ( 12 ) 

2 E g + , 2 E n + (2) 

V , X + (4) 
6 E g + , 6 S U + (6) 
2 n g , 2 n u (4) 

4 n g , X (s) 
X , 6 n n ( 12 ) 


Figure 3 is a plot of the potential curves of the 18 states of 0 2 that dissociate into 
oxygen atoms in the 3 P state. The discrete data points are shown along with the 
resulting spline fit. The corresponding curves for N-N, 0 + — O, and N + — N for dissocia- 
tion into ground state ion-atom pairs are shown in Figures 4, 5, and 6, respectively. 
Figure 7 shows the spline fits to the discrete data for 0-0 in their lowest states as well 
as the five excited states. (The discrete data points used to generate these plots have 
been omitted for clarity.) 




















Energy, an 




-0.30 -0.25 -0.20 -0.15 -0.10 -0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 


- 13 - 



Figure 5. 0 + (*S°)— 0( 3 P) Interaction Energies 
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Figure 7. Ground State and Excited State 0-0 Interaction Energies 
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IY. THE SKATR COMPUTER PROGRAM 


The SKATR program fits a spline curve to the discrete potential data and accepts any 
long range interactions of the general form 


V(R) = A^e A,R + 





For small values of R, the program determines the parameters of an exponential func- 
tion to agree with the spline Fit and the First derivative of the spline at the beginning 
of the discrete potential data. A typical input data set is shown in Table 2. 


Table 2. Potential Data for 0-0 State X 5 S g 

R (a 0 ) 

V(R) (a.u.) 

1.54700 

0.36800 

1.58700 

0.29400 

1.62700 

0.22100 

1.67700 

0.14700 

1.71700 

0.07400 

1.74700 

0.02825 

1.84740 

-0.05864 

1.87710 

-0.08165 

1.91560 

-0.10686 

1.96850 

-0.13417 

2.03840 

-0.16021 

2.13020 

-0.18084 

2.28195 

-0.19151 

2.47160 

-0.18084 

2.64110 

-0.16021 

2.81890 

-0.13417 

3.00070 

-0.10686 

3.17940 

-0.08165 

3.36410 

-0.05864 

3.40000 

-0.05456 

3.60000 

-0.03560 

3.80000 

-0.02141 

4.00000 

-0.01222 

4.50000 

-0.00333 

5.00000 

-0.00145 

5.50000 

-0.00078 

6.00000 

-0.00047 

6.50000 

-0.00030 

7.00000 

-0.00019 

FOR R<1.547, V(R) = 893.31 e (~ 5 0385R ) 

FOR R>7.000, V(R) = - 

-401.9 e (- 2 - 6002R ) 22 f 

R 6 
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The SKATR program uses the interaction potential to determine the effective poten- 
tial, 


V e (R) - V(R) + 


(/+y 2 ) z E 

k 2 R 2 




for increasing values of the angular momentum quantum number /. The maxima and 
minima of the effective potential are located and the turning points are determined for 
a large set of energies. This permits the phase shifts to be evaluated using the 
appropriate semiclassical formulas. The process is continued until the running sum 
for determining the collision cross sections satisfies a convergence criterion. 


Numerous tests were conducted to validate the code including those described below. 
Check on ^7j w kb5 Cohen [9] reports the results of accurate calculations of 7?j wkb /kr e 
for a variety of potential functions and various energies. The values reported by 
Cohen have a stated accuracy of five decimal places. For the Morse potential, 

with ar e = 3, D e = 0.324, and r e = 2.12, six cases were compared corresponding to 
E/D e =0.01 and 0.8 with /?/r e = 0.25, 1.00, 1.50. It was found that our results agreed 
precisely with those of Cohen to five decimal places. 

Check on Q n ; Equation (2) is applicable at high energies and the SKATR program 
can be applied at all energies (although it is quite inefficient to do so when the quan- 
tum effects are negligible). For the Morse potential (given above) and E>0.15 a.u., the 
energy is above the barrier of the effective potential (see Figure 2). Table 3 shows the 
results of the SKATR program compared with the classical results for the collision 
cross sections. (The mass of each collision partner was taken to be equal to that of 
nitrogen.) The values of the cross sections listed are obtained from Equation (2) 
using a code developed by Galant [10] and 8 is the difference between the classical and 
the quantum mechanical results. As the energy increases, the corrections become less 
significant as expected. 



Table 3. 

Comparison With Classical Calculation 

E (a.u.) 

Qi 


Q2 

^2 

Q 3 

<53 

mm 

59.67 

7.2 (-2) 

40.07 



life's® 

E9 

39.56 

3.8 (-4) 

38.12 



-gg'H 

1‘, 

17.41 

1.8 (-4) 

21.72 

[fifl 


* MK 


11.14 

-8.0 (-6) 

12.97 



HU 


Successful tests of Q n were also done using a Lennard-Jones potential of the form, 
V(r) =4 D e (r~ 12 — r -6 ) and comparing the results with O’Hara and Smith [11]. 
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V. RESULTS 

Principal results along with selected illustrative examples of intermediate results are 
presented in this section. A complete tabulation is contained in a report [36]. Using 
the input data set shown in Table 1, we evaluated the semiclassical equations to deter- 
mine the phase shifts and subsequently the collision cross sections Q lt Q 2j and Q 3 . 
Since the Q n are needed as functions of energy for the subsequent evaluation of the 
collision integrals (Eqn. (l)), this evaluation was done for a set of 40 E values from 
10"" 4 a.u. to 1.90 a.u. For each of the higher energy cases, the program typically 
required about 1500-2000 values of the angular momentum quantum number to 
assure convergence of the infinite summation of Equation (8) to five significant fig- 
ures. For very low energies, about 500 values of / were typically required for conver- 
gence. The running time on the Cray XMP-4/8 for the entire set of 40 energy values 
(corresponding to ~ 40,000 phase shifts r)[) was about two minutes for each state. A 
small subset of the results is shown in Table 4; the cross sections Q n are expressed in 
units of a 0 2 , where a 0 is the Bohr radius. 


TABLE 4. COLLISION CROSS SECTIONS FOR 0-0 STATE X 3 E “ 

E (a.u.) 

Qi 

Q2 

Q3 

1.50 (-0) 

5.00 (-1) 

1.00 (-1) 

5.00 (-2) 

1.00 (-2) 

5.00 (-3) 

1.00 (-3) 

6.521 (+0) 
1.291 (+1) 
4.319 (+1) 
4.680 (+1) 
5.642 (+1) 
7.856 ( + 1) 
1.349 (+2) 

6.470 (+0) 
1.477 (+1) 
2.980 (+1) 
3.433 (+1) 
4.721 (+1) 
5.854 (+1) 
1.026 (+2) 

8.860 (+0) 
2.012 (+1) 
4.748 (+1) 
5.226 (+1) 
6.629 (+1) 
9.104 (+1) 
1.591 (+2) 


In order to obtain the mean collision integrals for 0( 3 P) + 0( 3 P) collisions, each of 
the 18 states listed in the first column of Table 1 w r as processed as described above 
and combined according to their degree of degeneracy to yield weighted mean collision 
cross sections Q n . A sample of the results is shown in Table 5 and in graphical form 
in Figure 8. It was also of interest to investigate the contribution of each of the states 
to the mean collision cross sections. Table 6 shows the contribution of each 
0( 3 P)-0( 3 P) state to the mean cross section Q 1 . Similar tables were prepared for the 
mean values of Q 2) and Q 3 , and the process was repeated for all the excited states of 
Table 1. The results emphasize that the high spin, largely repulsive states are 
extremely important contributors to the collision cross-sections and hence to the tran- 
sport properties. 

For like ion-atom collisions such as N + — N, the phase shifts for each gerade and 
ungerade pair of states are combined to determine the cross sections according to 
Eqs.(4)-(7). Table 7 illustrates a sample of the resulting cross sections for an N + — N 
interaction. Note that Qj and Q 3 approach a value 2Q ex at high energies as expected; 
at low energies where the induced dipole polarization force becomes important, this 
simple relationship is no longer valid. 
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TABLE 5. MEAN 0( s P)-0( 3 P) CROSS SECTIONS 


E (a.u.) 



MEAN Qj 

MEAN Q 2 

MEAN Q 3 

8.3843 
1.6288 ( + 1) 
3.1451 (+1) 
3.9492 ( + 1) 
6.1801 (+1) 
7.5134 (+1) 
1.1891 (+2) 

9.8475 
1.5755 (+1) 
2.7872 (+1) 
3.3994 (+1) 
5.2321 (+1) 
6.1451 (+1) 
9.9196 (+1) 

1.2629 (+1) 
2.1214 (+1) 
3.8382 (+1) 
4.6898 (+1) 
7.2494 (+1) 
8.7187 (+1) 
1.4022 (+2) 


TABLE 6. CONTRIBUTION (%) OF EACH 0( a P)-0( a P) 
STATE TO WEIGHTED MEAN Q x 


ENERGY (a.u.) 


0.01 


0.005 


STATE 


1.5 

0.5 

0.1 

0.05 

1 

1 

2 

2 

1 

1 

1 

1 

1 

2 

2 

2 

3 

3 

5 

4 

3 

3 

2 

3 

3 

5 

5 

5 

7 

7 

6 

6 

9 

8 

8 

7 

10 

10 

8 

8 

1 

2 

3 

3 

2 

2 

2 

2 

5 

4 

6 

7 

3 

4 

5 

5 

14 

13 

10 

10 

15 

16 

14 

13 

2 

2 

4 

3 

6 

5 

5 

7 

15 

16 

14 

13 


0.001 


1 

1 

1 


b W 

c J E u - 

2 3 E ‘ + 

x 3 e- 

A 3 E u + 

2 3 E u + 

i 5 e u - 
1 % + 
2 5 E„ + 



C 3 A, 




TABLE 7. COLLISION CROSS SECTIONS 
FOR N + ( 3 P)-N( 4 S°) PAIR [ 4 E g + , 4 S U + ] 


E (a.u.) 


0.5 

3.2129 

(+2) 

1.1953 

(+1) 

3.2096 

(+2) 

1.5992 

(+2) 

0.15 

3.7586 

(+2) 

3.8853 

(+1) 

3.7606 

(+2) 

1.8772 

(+2) 

0.05 

4.2311 

(+2) 

6.3146 

(+1) 

4.2401 

(+2) 

2.0815 

(+2) 

0.005 

5.5623 

(+2) 

1.7158 

(+2) 

5.5643 

(+2) 

2.6724 

(+2) 

0.0005 

8.4220 

(+2) 

4.4297 

(+2) 

8.7678 

(+2) 

3.5775 

(+2) 
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The resonance charge exchange cross sections have been determined with the present 
code and the results are displayed in Figure 9 along with the asymptotic approxima- 
tion of Stallcop and Partridge [2] and available measured data from beam experiments 
[37,38]. Both the present "exact" results and prior approximate theoretical values 
agree with the measured data to within the experimental uncertainty. Hence, the 
present calculations give further support to the reliability of our calculated values of 
the long-range ion-atom interaction. Furthermore, a comparison of theoretical and 
experimental values provides a means to assess the contribution of excited state 
interactions. This contribution is expected to become appreciable at high energies; 
however, the results displayed in Figure 9 indicate it is small for the range of energies 
of the present work. 

Some principal results of the present study are summarized in Tables 8-11, which 
show the nine collision integrals (in units of square Angstroms^lO -16 cm 2 ) needed for 
the determination of the transport properties of a pure gas to second order or a gas 
mixture to First order. These tables represent the ground state O-O, N-N, 0 + — O, and 
N + — N collision integrals, respectively for temperatures ranging from 250 to 100,000 
K. Similar results for the excited state 0-0 collisions are contained in report [36]. 

The discrete data points were typically available only for a range of R from about 1.5 
to 10.0. Outside of this range, analytic extensions to the data were used. Experi- 
ments varying the short range extrapolation showed that the final results (for both 
the collision cross sections as well as the collision integrals) were relatively insensitive 
to the specific extrapolations for small R. Similar experiments for the long range 
extrapolation showed quite different results. Both the functional form as well as the 
particular parameters selected within a functional form had a significant effect on the 
collision cross sections for small values of the energy (e.g., energies below ~0.001 a.u.). 

Furthermore, as can be seen from Figures 4-7, several potential curves exhibit poten- 
tial barriers arising from avoided crossings of potential curves of other states with the 
same symmetry. For some states, the combination of this maximum and a well in the 
potential curve introduced a double barrier in the effective potential energy for low 
energy collisions. The resonance effects are expected to be small in this case for the 
scattering conditions of the present study. Consequently, a barrier treatment of the 
type described in Section II was not applied. The phase shifts were simply determined 
from the JWKB approximation. (This is equivalent to considering the region inside 
the outermost turning point to be inaccessible, as in a classical mechanics treatment.) 

In view of the above considerations, some of the cross sections at low energies may not 
be accurate because of uncertainties in the long range interaction energies for some 
states and our approximate treatment of the double barrier cases. However, the effect 
of these possible low energy cross section errors on the collision integrals is sharply 
diminished by the presence of the factor x 2s+3 in Equation (l). 

Similarly, the uncertainty in the tabulated values of the collision integrals due to 
extrapolation of the potential data at small R should not be significant for two rea- 
sons. In the first place, the scattering from the repulsive walls of inner low lying 




Collision Energy, eV 

Figure 9. Resonance Charge Exchange Cross Sections 
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TABLE 10, MEAN COLLISION INTEGRALS, cj 2 T^i,j), FOR 0 + ( 4 S°)-0( 3 P) COLLISIONS 

TEMP K II ct 2 Q( 1 , 1 ) |ct 2 Q( 1 , 2 ) | ct 2 Q[ 1 , 3 ) I a 2 0 ( 1 , 4 ) | ct 2 Q( 1 , 5 ) |<j 2 Q( 2 , 2 ) |a 2 Q(2,3) |a 2 n( 2 , 4 ) |ct 2 Q( 3 , 3 ) 
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states does not contribute to the collision integrals as heavily as the contributions 
from the high spin upper states. Secondly, although such extrapolations of the repul- 
sive wall may introduce uncertainties in the cross sections at high energies, the effect 
on the collision integrals is diminished due to the factor e~ x in Eqn. (l). 

Several collision integrals showing the interaction of oxygen atoms have been plotted 
in Figures 10, 11 and 12. The collision integrals a 2 fl(l,l) and a 2 fl( 2,2) for the 
interaction of oxygen atoms in their ground state from Table 8 are shown in Figure 
10. These curves are compared with our calculated results for interactions involving 
excited state atoms in Figures 11 and 12. Note that the results for the 0( 1 S)— 0( ! S) 
and 0( 1 D)— O^S) states clearly exhibit the high temperature trend expected of 
excited states. Since atoms in the excited state are "bigger", i.e., the outer electrons 
are further away from the nucleus, one expects that the interaction energy is larger 
and cr 2 fl would increase with excitation energy. However, Figures 11 and 12 show 
that the values of a 2 Cl for other excited state interactions do not all follow this trend. 
Part of this behavior is explained by the tendency of the repulsive wall of the poten- 
tial curve to be pushed to smaller R by curve crossing interactions. This results in 
smaller values of the collision integrals. Note also that the collision integrals of the 
0( l S)— 0( X S) interaction are small at low temperatures. This is expected because of 
the repulsive non-bonding nature of the interaction at large separation distances. 

We conclude our discussion of the results by displaying some plots of transport quan- 
tities that may be obtained readily from the tabulated values of the collision integrals. 
The viscosity for a gas of oxygen atoms in their ground states is shown in Figure 13. 
As can be seen by comparing Tables 8 and 9, the corresponding values for nitrogen 
are very close to these results. The thermal conductivity and self-diffusion coefficient 
of these two gases are compared in Figures 14 and 15, respectively. The binary diffu- 
sion coefficient for an N + — N mixture and for an 0 + — O mixture are also shown in Fig- 
ure 15. Note that the ion-atom curves lie considerably above the atom-atom curves; 
this is due to the fact that resonance charge exchange is an important mechanism for 
charge diffusion. The transport quantities shown in Figures 13-15 can be used as 
input for the determination of transport properties of mixtures such as that occuring 
in high temperature air. 


VI. CONCLUDING REMARKS 

We have developed a general method to permit a direct determination of transport 
collision cross sections from discrete potential data, rather than less reliable fits to 
analytical potential functions. Our semiclassical methods account for quantum 
mechanical effects such as resonance charge exchange, which is important at high tem- 
peratures, without the use of ad-hoc assumptions. The code can accomodate the case 
of multiple turning points; e.g., when the potential energy exhibits barriers which can 
arise from the non-crossing of potential curves for states of like symmetry (avoided 
crossings). 
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Figure 11. Collision Integrals a^^for 0-0 Interactions 







Thermal Conductivity, A, mcal/(cm-sec-K) 

• Figure 14. Thermal Conductivities for Oxygen and Nitrogen 
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Diffusion Coefficient, pD, cm' atm/sec 
Figure 15. Diffusion Coefficients for Oxygen and Nitrogen 



